perm filename FMSIG.F4[1,MUS]1 blob
sn#088194 filedate 1974-02-19 generic text, type C, neo UTF8
COMMENT ā VALID 00002 PAGES
C REC PAGE DESCRIPTION
C00001 00001
C00002 00002 DIMENSION W1(512),FUNC(512),SU(100),ST(550)
C00004 ENDMK
Cā;
DIMENSION W1(512),FUNC(512),SU(100),ST(550)
GO TO 60
40 CALL HYDPOG(2)
60 TYPE 50
50 FORMAT(' TYPE A,CF,MF,I1,I2,STEPS'/)
ACCEPT 100,A,FC,FM,XI1,XI2,STEPS
100 FORMAT (6F)
CALL TYPLOC(-300,-512)
CALL DPYSET(1,SU,100)
CALL DPYBRT(1)
CALL AIVECT(0,0)
CALL ALINE(0,-12,0,12)
CALL ALINE(256,-12,256,12)
CALL ALINE(-268,128,-244,128)
CALL ALINE(-268,256,-244,256)
CALL ALINE(-268,-128,-244,-128)
CALL ALINE(-268,-254,-244,-254)
CALL ALINE(-264,0,256,0)
CALL ALINE(-256,-256,-256,256)
CALL DPYOUT(1)
XI=XI1
XI1=(XI2-XI1)/STEPS
N=STEPS+1
CALL DPYSET(2,ST,550)
CALL DPYBRT(2)
Y=6.2832/512.
DO 400 N1=1,N
CALL CLRPOG(2)
I=-1
DO 200 J=1,512
X=J-1
FUNC(J)=A*(SIN((FC*X)*Y-XI*(COS((FM*X)*Y))+XI))
IY2=FUNC(J)*256.
IF(I.EQ.0)GO TO 150
CALL AIVECT(-256,IY2)
I=0
GO TO 160
150 CALL SVECT(1,IY2-IY)
160 IY=IY2
200 CONTINUE
CALL DPYOUT(2)
XI=XI+XI1
400 CONTINUE
TYPE 403
403 FORMAT(' TYPE 0 TO CONTINUE'/)
ACCEPT 300,M
300 FORMAT (I)
IF(I.EQ.0)GO TO 40
END